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Abstract 

We analytically study phonon transmission and localization in random superlat- 
tices by using a Green's function approach. We derive expressions for the average 
transmission rate and localization length, or Lyapunov exponent, in terms of the su- 
perlattice structure factor. This is done by considering the backscattering of phonons, 
due to the complex mass density fluctuations, which incorporates all of the forward 
scattering processes. These analytical results are applied to two types of random 
superlattices and compared with numerical simulations based on the transfer matrix 
method. Our analytical results show excellent agreement with the numerical data. A 
universal relation for the transmission fluctuations versus the average transmission is 
derived explicitly, and independently confirmed by numerical simulations. The tran- 
sient of the distribution of transmission to the log-normal distribution for the localized 
phonons is also studied. 

PACS numbers: 62.65. +k, 63.20.-e, 68.65.+g 
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I. INTRODUCTION 



The long-standing interest in Anderson localization has greatly increased following the un- 
derstanding that it is common to both electrons and classical waves in disordered mediaS. 
For the localization of acoustic phonons, an experimentally relevant quantity to study would 
be the transmission rate. The well-known basic feature is that phonons will not propagate 
through a medium with a large amount of randomness. In particular, in one dimension any 
disorder is strong enough to induce the exponential localization of eigenmodes, i.e., their 
amplitudes decay exponentially along the medium with a rate of decay called the localiza- 
tion length. So far, phonon transmission rates have been measured in synthetic multilay- 
ered systems, or superlattices (SL's), with both periodic and quasiperiodic order!. Similar 
experiments with randomly-layered systems would provide important information on the 
localization properties of high-frequency acoustic phonons. 

In a previous paper!, we studied phonon propagation through random SL's by means 
of the transfer matrix method. The random SL's considered are the multilayered systems 
where two kinds of basic blocks of materials (which may or may not have internal structure) 
are stacked at random. Each sample has its own realization of randomness in the order 
of constituent layers, and the phonon transmission rate versus frequency shows a fine spiky 
structure specific to the particular realization of disorder present in that given sample. This 
fluctuating transmission rate can be considered to be a fingerprint of the sample, like the re- 
producible conductance fluctuations found in electronic mesoscopic transport. The ensemble 
average of the transmission rate over possible configurations of the constituent layers smears 
out the fine structure, but still leaves global features of transmission dips and peaks. The 
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dips arise from phonon localization due to the interference among backscattered phononso. 
The peaks arise from the resonance occurring when appropriate matching conditions are 
satisfied between the phonon wavelength and the thickness of a basic layer. 

A significant result of Ref. 3 is that there exists a remarkable correlation between the 
ensemble- averaged reflection rate and squared SL-structure factor calculated analytically. 
The latter describes the sum of the phonon amplitudes reflected from SL interfaces. The 
smaller contributions coming from multiple phonon reflections have been neglected. We find 
that at the frequencies where the maxima (minima) of the structure factor are attained, the 
reflection rate exhibits peaks (dips). 

The purpose of the present study is to establish quantitatively the relation between the 
phonon transmission rate and the SL-structure factor. Also, we examine in some detail the 
localization characteristic of phonons injected into random SL's. Our study is based on the 
Green's function method which has originally been applied to the electronic conductivity in 
one-dimensional metalsS. 

In section II, we model the random SL's and formulate the phonon transmission rate 
in terms of Green's functions. In section III, we derive integral equations for the Green's 
functions by introducing complex mass-density fluctuations which induce the backscattering 
of phonons. The solutions of these integral equations, satisfying appropriate boundary 
conditions, are given in terms of the S-matrix associated with the problem. In section IV, 
we show that the S-matrix elements can be written in terms of the SL-structure factor, which 
describes the phonon amplitude backscattered from the complex mass-density fluctuations. 
The Born approximation is employed for the scattering and the explicit expression is given 
for the phonon mean-free-path. In section V, we derive an analytical expression for the 
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average transmission rate (T), as a function of the system size divided by the phonon mean- 
free-path. We also apply this formula to two types of random SL's and compare the 
analytical results with the numerical simulations based on the transfer matrix method. In 
section VI, the Lyapunov exponent, or the inverse localization length, is obtained by using 
a differential recursion relation. The Lyapunov exponent is given by the reciprocal of 
the mean-free-path and exhibits a good agreement with the numerical simulations. The 
universality of the transmission fluctuation is discussed in section VII. In section VIII, we 
study the probability distribution of the transmission rate and confirm that the log-normal 
distribution is applicable to the localized phonons. A summary and discussions are presented 
in the last section. 

II. Model 

In the present study, we consider the phonon transmission through two types of random SL's 
with different structures in their constituent layers. The first class of SL's, which we refer 
to as single-layer SL's, consist of two kinds of unit blocks as illustrated in Fig. 1(a); that is, 
each unit block is composed of a single material (A or B) with a definite thickness (dA or 
d B ). The second class of SL's, denoted as double-layer SL's, also consist of two kinds of 
unit blocks but each one has an internal structure (or basis) as shown in Fig. 1(b). The 
unit blocks of the latter class (specified by indices 1 and 2) are composed of two different 
kinds of layers of materials A and B with definite thicknesses (d± t A and d 1B in the first block 
and G? 2 ,A and g^.b in the second block). We assume that the order of the materials A and 
B is the same in these building blocks of double-layer SL's. Single-layer random SL's are 
constructed by stacking these two kinds of single-layer unit blocks at random with equal 
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occurrence probabilities, i.e., 0.5. A random stacking of the two types of double-layer unit- 
blocks creates a double-layer random SL. It should be noted that phonons cannot recognize 
the existence of the interfaces between consecutive AA and BB layers which occur very often 
in the single-layer SL's. So, it is convenient to define a "segment" which stands for any part 
of a single-layer SL consisting of a consecutive series of the same kind of blocks. In short, 
a segment is a part of a SL consisting of a single material A or B. In any double-layer SL, 
however, different materials always meet at any interface. Hence, for the double-layer SL the 
segment is only a part of a unit block. Hereafter, we call a segment consisting of material 
A (B) simply as an A (a B) segment. 

If the SL's considered are composed of two kinds of materials, A and B, they can be 
modelled as an alternating sequence of segments of materials A and B as described above. 
The thicknesses of the segments, denoted by D A and D B , are distributed according to the 
prescribed structures of the basic building blocks and the probability for the successive 
stacking of the same kind of material. For simplicity, A and B are assumed to have the 
same stiffness constant \i and are distinguished by their mass densities p A and p B . We also 
assume that a random SL is sandwiched between a substrate and a detector composed of 
material A. 

Now, we consider the situation where an acoustic wave excited in the substrate propagates 
through a random SL and is observed at the detector on top of the SL. We assume the 
propagation to be perpendicular to the interfaces of the building blocks. No mode conversion 
among phonons at the interfaces is considered. This is valid if the interfaces have a mirror 
symmetry. As the randomness is due to the stacking order of blocks and is present only along 
the propagating direction (x-direction), the system we discuss is essentially one-dimensional. 
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The basic wave equation governing the motion of the displacement u(x, u) at an angular 
frequency uj is 



where the mass density p(x) takes on the value of either p& or p B depending on the position 
x (0 < x < L) in the SL of length L. Also, p(x) is taken to be pA in the substrate (x < 0) 
and detector (x > L). We define 5p = ps — Pa- For convenience, we assume pa < Pb but 
the results are also valid for pa > Pb- 

We formulate the transmission in terms of the retarded and advanced Green's functions 
associated with Eq.(|l]): 



where uj± = uj ± id and 5 is an infinitesimal positive number. We define the transmission 
rate by the ratio of the acoustic Poynting vector in a random SL to that in a homogeneous 
system with the mass density pa- The unperturbed Green's functions satisfy the equation: 



(p(x)uj 2 + pdl)u(x,uj) = 0, 



(1) 



(p(x)a;| + iidi£)Gu±(x, x') = 5(x - x'), 



(2) 



{p A ujl + pdl)Gl ± {x,x')=5{x-x') 



(3) 



and the solutions are 



G° w± (x,x') = t 



±iu}\x—x'\/cA 



(4) 



2p A c A uj 



where ca is the sound velocity in the A material which is given by 




(5) 



and cb is the sound velocity in the B material, 




(6) 



6 



Thus, the transmission rate T(L, uj) is written with G u ± as 

T(L,cu) = (2p A c A tu) 2 G ul+ (L,0) G u -(L,0). 



(7) 



III. One-particle Green's function 

The mass density fluctuation 5p(x) = p(x) — p A present in the SL causes the forward and 
backward scatterings of phonons incident upon the SL from the homogeneous substrate. A 
key observation is that the effect of the forward scattering processes can be exactly summed 
up. The Green's function G° incorporating all of the forward scattering processes is related 
to the unperturbed Green's function G° as 

G° w+ (x,x') = G w+ (x,x') e i«<ff«<*-*')[*(*) - ( 8 ) 

where the phase $ is defined by 

= ir— r^p(y)dy. (9) 

2p A c A Jo 

Now, considering only backscattering processes, we can express the solution of Eq. (^|) in terms 
of the phase-modulated Green's functions, G° + 's. For an even number of backscattering 
processes, and for x > x', the retarded Green's function satisfies 

G u+ (x,x') = Gl + {x,x') G^ + (x,x 1 )(-5p(x l )u 2 )Gl + (x l ,x 2 )(-8p(x 2 )u 2 ) 
xGuj + (x2, x')9(x2 — xi)9(x — Xi)9(x2 — x')dx\dx2. (10) 

In order to express Eq.flKp with the bare Green's function + , we multiply both sides of 
Eq.(|lO"D by exp[— z($(x) — and introduce 

G u+ (x,x') = G^(x,x')e- 1 ^-^. (11) 



Thus, we obtain 



GW(x,x') = Gl + (x,x') +j o j Q G° UJ+ (x,x 1 )(-ax l )u J 2 )G UJ+ (x l ,x 2 )(-C(x2p 2 ) 

xG UJ +(x2, x')9(x2 — xi)9(x — Xi)9(x2 — x r )dx\dx2 1 (12) 

where ((x) denotes the complex mass density fluctuation defined by 

C(x) = 6p(x)e- 2i ^ x \ (13) 

As the complex mass density fluctuations, (* and (, transfer a phonon with a wave number 
k(= uj/ ca) to — k and vice versa, it is convenient to adopt a 2 x 2 matrix representation for the 
Green's functions, > where the indices i, j — 1 and 2 correspond to the phonon states 

with k and — k, respectively. Defining the diagonal matrix elements for the unperturbed 
Green's function, g® + , as 

\9°+{x,x')] n = G^ + (x, x') 9(x - x') 

= 9(x - a; ') e M*-*')/<* n 4 ) 

and 

[gl + U{x,x') = Gl + (x,x') 9{x' - x), (15) 
we re-express Eq.(|T2|) in terms of the matrix elements of g w+ and g^ + as 
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[g u+ (x,x')] n = [#° + (x,x')]ii + j j [^° + (x,Xi)] n (-C(xi)w 2 ) [^ + (xi,x 2 )] 2 

x (-C(x2)^ 2 ) [9w+(x2, x')} 11 dx 1 dx 2} (16) 

where is identical to G w+ defined by Eq.(|TT]) and the indices indicate that it begins 

and terminates with [g° + ]n. The other diagonal element [g^+]22 is obtained from Eq. (|I6"D by 
interchanging ( and (* and the indices 1 and 2. 



We can treat the case of odd numbers of backscattering events in the same way. The 
off-diagonal elements [5^+] 21 defined by 

[g^+(x, x')\ 2 x = G u+ (x, x') exp[ i ($(x) + $(x'))], (17) 

can be expanded in terms of the complex mass density fluctuations ( and (* as follows: 

[gu+{x : x')\ 21 = / [gl + {x, xi)} 2 2 (-C*(xi)^ 2 ) [g° + (x u x'^udxx 



+ 



L r L 



JO 
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x(-C*(x 1 )u ) [g^x^x'^ndXi dx 2 . 



The expression for the element [<7uH-]i2 defined by 



[gcu+(x,x% 2 = G u+ (x,x!)exp[- i ($(x) + $(x'))] 



(19) 



is obtained from Eq.(|18|) by interchanging the ( and (*, and also the indices 1 and 2. Now, 
we can rearrange the equations for g^+, by employing Pauli's spin matrices 03 and = 
<j\ ± 102 , to give 



2p A oo [tc A a 3 d x + uj + }~y \<T W £(x) + cr^C(x)} \ ^+(x, x') = 8(x - x'). (20) 



To solve Eq.(^), we follow the mathematical method developed by Abrikosov and Ryzhkin 
for treating the electronic conductivity problem in one-dimension^. 

First of all, we introduce the 2x2 S'-matrix which obeys the following equation : 



2p A u [ic A cr 3 d x +u]-^r <t (+) C(i) + a [ - } e(x) \ S u (x, x') = 



Ha*/ 



(21) 



As the terms with complex mass density fluctuations ( and (* do not commute with each 
other, the S'-matrix can be expressed with the position-ordering operator T x as 

1 



S w (x,x') =T x exp\i— / {cr 3 — - — 
I c A Jx> 1 Ap A 



<r (+) C(ari) -<T (_) C*(a;i)]}tiPi 



(22) 



The operator T x rearranges a product of position-dependent operators so that the operators 
must be placed from right to left in the order of increasing x. For x > x\ > x', the retarded 
Green's function g u+ (x^) is related to g u+ (xx,x r ) by the expression 

gu>+{x,x') = S^x.x^g^x^x'). (23) 

The above relation also holds for x' > x > x\. However, for x > x' > x±, the relation between 
gv+(x, x') and gu+(xi, x') is not given by Eq . (p3j) because g u +(x, x') has an additional term 
arising from the ^-function in Eq.(^). After a bit of algebra, we obtain 

i 

g u+ {x, x') = S^ix, x^g^xx.x') - - S u (x, x')a 3 . (24) 

2p A c A u 

Equations (|20|) , (|23"D and fl24]) hold for the advanced Green's function as well as for the 
retarded Green's function. 

Now, we impose the boundary conditions satisfied by the retarded Green's function on 
Eqs.fl2"3]) and fl2~j|). Noting that the perturbation series of Eq.(|TjJ) for [g w+ {x, x')] u begin and 
terminate with [<7° + ]n and that all terms in the iterative solution [g® , (x, x')]n vanish in the 
limit of x' — > oo or x — > — oo, the boundary condition for \g u+ (x,x')]n becomes 

[9u+(x, x ')]n — ¥ at x — > — oo or x' — > oo. (25) 

In the same way, we obtain the following boundary conditions for the other elements of the 
retarded Green's function: 

[guj+{x,x')]22 ^0 at x — » oo or x' — > — oo, (26) 
[fiw( x > x ')]i2 — ► at x — > — oo or x' — > — oo, (27) 
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\g,j+(x, x')]2i — > at x — > oo or x' — >■ oo. (28) 
Taking the limit of £1 — ► —00 in Eq. (|24]) , we obtain from the boundary conditions (|25|) to 



\g u+ (x, x')] a /3 = [S u (x, -oo)] q2 [^+(-oo, x')] 2 (3 - [S w (x, x')a 3 } a3 9(x - x'). (29) 

2p A c A u 

If we take the limit of x — > 00 and put a = 2, the left hand side of Eq.(p9[) vanishes. Thus, 
we find 

i 1 

[^+(-00, x%/3 = - TF - rr- [5 w (oo, x')a 3 } 20 . (30) 

2p A c A u [bv{oo, -00JJ22 p 

Substituting Eq . (|30"D into Eq.(p9|), we can express the retarded Green's function as follows: 

K+(x, x ) a/3 = 4 r— ; r= 5 w (x, x )a 3 \ a g 9{x - x ) \ . (31) 

2p j4 c j4 w [ [60,(00, -00 J J 22 M J 

Similarly, the advanced Green's function can be expressed as 

^ w -(x,x) a/? = <^ — r= [^(x,x)cr3j a/3 ^(x-x) f . (32) 

2p A c A uj { [6^(00, -oojjn J 



Equations ( pT|) and (p2|) are used to derive the transmission rate. 

To proceed further, we introduce an interaction representation of the S-matrix with 
respect to uj 

S w {x, x') = e luJxa3/cA S{x, x')e- WfT3/cA . (33) 

Explicitly, 



S(x,x') = T.exp j~ [aM&M ~ * H C(si)] dxA , (34) 



X 



where 



= ia( + )C(x)e- 2 ^ x/c ^. (35) 
11 



It should be noted here that the S-matrix S(oo, — oo) in the interaction representation is 
identical to S(L, 0), or 

5(00,-00) = S(L,0), (36) 

since ((x) and (*(x) vanish outside (i.e., x > L and x < 0) the SL. 

Next, we divide the system into a set of slabs with thickness A and perform the integral 
of Eq. fl34|) . To do this, we define Kj by 



—1- 



2p A c A 



CM d y> 



(37) 



where Xj (= jA) denotes the position of interface between the jth and (j + l)th slabs. In 
terms of Kj and its complex conjugate the integral in Eq.([!4]) is replaced with a sum and the 
S-matrix becomes a product of the factors containing Kj and k* at different slabs. Keeping 



terms to second order in Kj and k*, we obtain 



rr/A 

n 

j=x'/A 



I 



\ 



1 + 2 1 K i 1 2 



l + ik,l 2 



(3? 



2l^l J 



We note here that the determinant of the S'-matrix is unity to order |k| 2 . 



IV. Backscattering Rate and SL Structure Factor 

The quantities {Kj} constituting the S-matrix are a set of fundamental random variables 
which vary from sample to sample. The statistical properties of {Kj} are expected to reflect 
the structures of the basic building blocks of random SL's. Equation ( |37"|) can be rewritten in 
terms of the mass density fluctuation Sp(x) and sound velocities in the constituent materials 
A and B. Substituting Eqs. (|13D and (^) into (|37D , we have 

12 




The sound velocity c(x) has value c A in A and cb in B. In deriving Eq. (|39|) , we have used 



c B ~c A [l + 



2p A 



(40) 



Next we define the average of variables {/%} in a random SL by 



A 



L/A 



T £ K r 



(41) 



3=0 



Putting Eq.([39|) into Q41D , we can relate Kj to the structure factor Ssl(L,lv) of a random SL 
as follows: 



A 



o c(j/' 



where 



A6j)Cb_ 
L Ap A c A 



S S l(L,u) 



(42) 



S S l{L,u) = (— ^] / dy5p(y)exp 



y 2u 



o c(y' 



(43) 



By performing the integral for a random SL which begins with segment B, the structure 
factor Ssl(L, uj) is written as the superposition of amplitudes of phonons reflected once from 
each interface 

M [ 2j-2 

S SL (L,co) = pl-e-^exp -^^ , (44) 

jf=l m=0 

where m = 2uD m /c m (put O = 0), -D m is the thickness of the mth segment, c m is the sound 
velocity, and M is the number of B segments. The expression of the structure factor for a 
random SL which begins with segment A is similarly given by 



M 



S S l(L,u) = £(l-e 



-i8 2j 



cxp 



2j-l 
m=0 



(45) 
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The size of the random SL, L, is the sum of the thicknesses of its constituent segments 

N 

L = £ An, (46) 

m=l 

where N is the total number of segments in the SL. 

Now we consider the ensemble average (Tcj) of (oc Ssl) over possible realizations of 
random sequences of the segments. However, (Kj) vanishes because Ssl(L,u) is a complex- 
valued random variable. The ensemble average of the product KjK% also vanishes for j ^ k 



by the same reason. The non- vanishing term (|«j| 2 ) can be given in terms of the structure 
factor Ssl(L,u) as follows: 

= j# (\S 8L (L,u,)\*). (47) 

Here we note that, to lowest order in the mass density difference dp, the prefactor Spcs/^PA ca 
in Eq.(|4*^) is equivalent to the reflection coefficient Rba of an acoustic wave incident on the 
interface from the material B to the material A, 
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Z B + Z a 

= -R AB = R. (48) 
Here, Za and Z B are the acoustic impedances and are defined by Za = PaCa and Z B = PbCb- 



The quantity (|/%-| 2 ) is related to the phonon backscattering rate T caused by the complex 



mass density fluctuation ((x). This can be seen by expressing (|k., | 2 ) in terms of ((x). From 



Eqs.(J55|) and Q57D, (|^j| 2 ) can be rewritten as 

rL 2 \ 

(49) 



(|«,f) = A 



a; 2 



Ap A c A L 



o 



L C(s)e- 2i ^ /c - 4 da; 
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The backscattering rate T is calculated in the Born approximation as 



r(L,u) 



^p\c A L \ 



C(x)e~ 2iuJx/cA dx 



(50) 



Comparing Eqs. fl47|) , (f49|) and (|50|), the mean-free-path £ limited by backscattering is 

<(L;1J) - 1 = E(^ = ^<I^(^)I 2 >, (51) 



ca 



so we can express (|/%| 2 ) in terms of £ as 



<hf> 



(52) 



£(L,^)' 

In the rest of this section, we analytically calculate (\Ssl{L,uj)\ 2 ) for the random SL's. 
For an odd N, the first and last segments of a random SL are the same. For an even N, 
a random SL which begins with an A (B) segment terminates with a B (A) segment. For 
a double-layer SL, N is always even and it begins with an A segment and terminates with 
a B segment. If an end segment of a SL is an A segment, there is no phonon reflection at 
the boundary between the random SL and the adjacent substrate or detector (assumed to 
be made of A material). In this case, we regard the end segment as a part of the substrate 
or detector. Hence, we may consider the random SL's whose end segments are B and take 
N = 2M — 1. Thus the length of the random SL becomes L = ((D A ) + (D B ))M for large M. 
Here (Da) and (Ab) denote the averaged thicknesses of the A and B segments, respectively. 

By using Eq. (fP|) , the squared SL structure factor becomes 



\S SL (L,w) 



Wi_ e -<fe»-i) 

71=1 

M m-1 



m=2 n=l 



2m-2 

-i £ Q 3 

j=2n 



+ c.c. (53) 



Since there is no correlation between the lengths of adjacent segments, (9,, 's associated with an 
A segment (9j = 9 A = 2u D A / ca) and B segment (8j = Ob = 2uDb/ cb) are independent 
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variables. From these considerations, the ensemble average of the squared SL structure factor 
per segment, I s , defined by I s = (\Ssl(L, cu)\ 2 )/N, becomes! 



I s (L,u) 



Re 



;i-e A )(l-e B ) 
1 - e A e B 



{Da) + {D B ) 



Re 



l-e B \ 2 
1 - e A e B / 



(54) 



where €a = {e~ WA ) and e B = {e~ tdB ). Equation (|54"D gives the backscattering rate in the 
random SL's (through Eq.(|5"l~D). The second term of Eq.(^) depends on the system size L 
and vanishes in the limit L — > oo. 

Next we consider the case for €a = 1 and eb 1. This is one of the resonance conditions 
under which phonons can be transmitted perfectly through the SL. This condition is satisfied, 
for instance, when twice the thickness of A layer or A block is an integer-multiple of the 
wavelength of phonons. In this case, the SL can be effectively regarded as a bulk layer of 
material B. However, the phonon reflections at both ends of the SL still exist. Under this 
condition, only the second term of Eq. ([54]) survives. Thus, we can attribute the second term 
to the scattering at the SL-substrate and the SL-detector boundaries, which is important at 
frequencies close to the resonance. However, for 7^ 1 and e# = 1, i.e., at the resonance 
for B segments, both terms in Eq.(|54]) vanish. This asymmetry originates from the fact that 
the SL is sandwiched between materials of type A. 

Based on the above considerations, we may separate the mean-free-path i into two parts: 
one limited by scattering within the SL, £sl, and the other limited by boundary scattering, 
4- Thus, 

(55) 



— f- 1 4- f- 1 



where 



2R 2 



{D a ) + {D B ) 



Re 



;i-e A )(l-e B ) 
1 - e^es 



(56) 
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and 



£ b {L,uj) 



2R 2 



Re 



1 - e^es 



(57) 



V. TRANSMISSION RATE 



Noting that g w ±'s differ from G^^'s by phase factors only [see Eqs. ([□]), ([16]) , ([TT| ) and 
(|19|)1, we can re-express the ensemble averaged transmission rate (T) in terms of the Green's 
functions gu±(x,x') as 



(T(L,u)) ={2 PA c A uf{{ gul+ {LMii [0 U -OM)]22> 



(58) 



Substituting Eqs.([H]), (|32|) and (|36"D into (p8|), the average transmission rate can be written 



as 



(T{L,u>)) 



(59) 



(60) 



\[5(L,0)] n [S(L,0)] 22 /- 
For convenience, we introduce a set of functions {A n (L)} which are defined by 

A M I {[S{xMi2 [S(x,0)]2i} n \ 
n{ ' " \ {[£(*, 0)] n [S(x,0)]22} n+1 / ' 

The averaged transmission rate (T(L,uj)) is equal to A (L). From Eq. (|38|) , we obtain the 

relation between S(x, 0) and S(x — A, 0) 

/ . . \ 



S(x,0) 



1 1^12 



1 + 



V 



K 



l + |l«| 2 



S{x- A,0), 



(61) 



where k is defined by Eq. fl39|) with the integral over the interval [x — A, a;]. We substitute 
Eq . (|B~TT) into ( |60"1) and expand A n (x) in a power series in n keeping terms up to |k| 2 . Thus 
we get a difference equation for A n (x) 



A n (x) - A n (x - A) = (\k\ 2 ) 
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{n 2 A n ^(x - A)+(n + l) 2 A n+1 (x - A) - [n 2 + (n + l) 2 ] A n (x - A)} . (62) 



In Eq.(|62|) we have replaced the average (|k| 2 ) with (|k| 2 ). The average (\k\ 2 ) coincides 



with (\nj\ 2 ) of Eq. ([52]) . Hence, substituting Eq.(|52"D into (|62|) and taking the limit of A — > 0, 
we obtain the following equation for A n (t) : 



dAJf) 
dt 



n 2 A n _i(t) + (n + l) 2 A n+1 (t) 



n 2 + (n + l) 2 



A„(t), 



(63) 



where we have set x = L and changed the variable to t defined by t = L/£. Next, constructing 
a generating function A(z, t) given by 



A(z,t) = £ A n (t)z n , 

n=0 

we can derive the partial differential equation satisfied by A(z,t), 



(64) 



dA 
~dt 



z(l-z) 



o 2 a 



dA 



gz2 + <1-,)<1 -3,)- - (1-,)A 



(65) 



At £ = or L = 0, A„ = S n fi, since the S'-matrix becomes a unit matrix. Hence, the 
boundary condition for the generating function is A(z,t = 0) = 1. The solution of Eq . (|65|) 
satisfying this boundary condition isi 

s f°° 27rAtanh7rA 1 /l 1 , z \ (i ,±2\ f „ . . 

jo cosh7rA l — z \2 2 1 — zJ 

where F(a, 0,j;y) is the associated hypergeometric function. Since the ensemble averaged 
transmission rate (T(L,u)) is given by A(z = 0,t), the final expression for (T(L,u)) yields 



{T{L,u>)) 



00 2 7rAtanh7rA 



coshvrA 



cxp 



- - + A 2 t 



dX. 



(67) 



It should be noted that (T) is a monotonically decreasing function of t. We also note that 
(T) depends on phonon frequency only through the phonon mean-free-path I. 



A. Single-Layer Random SL's 



We shall now apply the analytical result (|67| ) for (T) to real systems consisting of AlAs- 
GaAs multilayers. First, we consider single-layer SL's constructed by randomly stacking 
two kinds of basic blocks, one being an AlAs (A) single layer with thickness cIaias and the 
other a GaAs (B) single layer with thickness dc a As [Fig. 1(a)]. The interfaces are assumed to 
be (lOO)-plane. The stiffness constant p of our model is the average of the C44 moduli for 
AlAs and GaAs. This approximation is good since C44 (= 5.89 x 10 11 dyncm~ 2 ) for AlAs 
is very close to the value for GaAs (5.94 x 10 11 dyncm~ 2 ), but paias (= 3.76 gcm~ 3 ) and 
PGaAs (= 5.89 gcm~ 3 ) are quite different. 

The ensemble average of the phase factors, €a = €aia s and e# = ^GaAs, are calculated as 

€AlAs = 2-e~ ia) €GaAs = 2-e~ ib) ^ 

where a = 2udAiAs/cAiAs, and b = 2ud Ga A s /c Ga A s , respectively. The average thicknesses of 
the segments are also evaluated as 

(Daias) = 2dAiAs, (DcaAs) = 2d GaAs . (69) 



Substituting Eqs.(p8|) and ( p9|) into (|5q) and (pf\ ), we obtain the expression for the mean- 



free-path £. The dominant part, £sl, of £ is 



-1 



AR 2 (1 - cosa)(l - cos 6) 



(70) 



SL (Daias) + {D GaA s) 3 - 2 cos a - 2 cos b + cos(a - b) ' 
where R = Raias,GclAs = {Zaias - Z GaAs ) / {Z A iAs + Z GaAs ) = 0.091 for a transverse mode. 

To verify the analytical results we have derived, numerical simulations based on the 
transfer matrix method are also made. Figure 2(a) compares (T) versus frequency up to 
1.3 THz for transverse phonons propagating along the growth direction of random SL's with 
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1000 blocks. The thicknesses of the blocks are chosen as (Iaias — dc a As = 34A, which 
correspond to twelve mono-atomic layers. The solid line indicates the analytical result. The 
open circles are the numerical data obtained by averaging over 100 random SL's with different 
sequences of blocks. The transmission rate of the periodic SL consisting of alternating AlAs 
and GaAs blocks is also shown with a dotted line for comparison. 

The average transmission (T) shows alternating prominent dips and sharp peaks. The 
large dips in transmission reflect the localization of phonons due to interference effects among 
the scattered waves. The transmission is perfect or almost perfect at the frequencies satis- 
fying the resonance conditions (£sl = °o). As the sound velocities of the TA-mode waves 
are cc a As = 3.32 x 10 5 cm/s and caias = 3.97 x 10 5 cm/s, the resonance frequencies are 
v GaAsn = n x 0.489 THz for a GaAs single-layer and Vaiasu = n x 0.583 THz for an AlAs 
single-layer. 

We see in Fig.2(a) the slight reductions of (T) from unity at v — v MAs n (n = 1, 2). These 
are due to the scattering of phonons at both ends of random SL's. Although i^ 1 goes to as 
L — ► oo, the deviation of (T) from unity survives even for infinitely large systems, because 
L/£b is a constant regardless of the system size L. 

Finally, we compare (T) with the behavior of the phonon transmission in the correspond- 
ing periodic SL. We see that the minima of (T) always occur in the band gaps of the periodic 
SL. From Eq.([7(]) we can show that the frequencies at the minima of (T) coincide with the 
Bragg frequencies! of the corresponding periodic SL satisfying sin [(a + b)/2] = 0. This 
apparent similarity in the transmission dips between the random and periodic SL's is quite 
interesting. Of course, the periodic and random SL-systems are very different. The periodic 
SL has a finite band width. However, as will be discussed below, the width of the pass bands 
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in the random SL's where the localization length is longer than the system size tends to zero 
for L — > oo. 



B. Double-Layer Random SL's 

The second type of random SL's to be studied here are the double-layer random SL's whose 
basic building blocks are shown in Fig. 1(b). In contrast to the single-layer random SL's, 
a characteristic feature of these SL's is that the interfaces between adjacent blocks never 
disappear even if we stack two kinds of basic blocks at random. Denoting the thicknesses of 
the constituent layers in the first blocks by d^AiAs and di tG aAs, and those of the second block 
by d2,AiAs and d 2tGa As, the ensemble- averaged phase factors caia s and e Ga As are given by 

g— iai _|_ g— io,2 



tAlAs — 



e -ibi + e -ib 2 



tGaAs = " , (71) 

where a 1{2 ) = 2ud 1 ^ 2 ),AiAs/cAiAs and &i (2 ) = %ujd 1{2) ^ GaAs / c GaAs . The average thicknesses of 
the segments are 

/ n v di AlAs + d 2 ,AlAs i \ d\ t GaAs + d 2 ,GaAs ,„n\ 

{■L>AlAs) = ^ ' K^GaAs) = ~ • (72) 

Hereafter we consider the case where d^AiAs = d 2j AiAs = dAiAs for simplicity. For this type 
of random SL, becomes 



! R 2 (1 — cosa)(l — cos(6 — c)) 

' SL ~ (Daias) + {D G aAs) 3-2 cos(a + b) - 2 cos(a + c) + cos(6 - c) ' 



(73) 



where a = a\ = a 2 , b — bi and c = b 2 . 

The evaluation of (T) for double-layer random SL's is performed by assuming 100 blocks 
(M = 100) with d A i As (= d lt AUs = d 2)A iAs) = 17A and d hGaAs = 20A and d 2)GaAs = 42A. 
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These pairs of building blocks are the same as those used by Merlin et al@ for the fabrication 
of quasi-periodic SL's. 

Figure 2(b) shows the analytically calculated (T) together with the numerical data. 
Again, the agreement between these two calculations is excellent. The frequency depen- 
dence of (T) looks quite different from Fig. 2(a) for the single-layer SL's due to the difference 
in the structures of the building blocks. In the present case, as can be seen from Eq.fl73|). 
the resonances occur for a = 2mr for the AlAs layers and b — c = 2mr for the GaAs lay- 
ers. The first condition implies that this resonance happens when an integer-multiple of the 
wavelength becomes equal to twice the thickness (Iaias of an AlAs layer. Numerically, the 
corresponding resonance frequency is given by Vaiasu = nxl.17 THz. The second condition 
means that the thicknesses of GaAs layers in the first and second blocks become effectively 
identical at the frequencies satisfying this condition (numerically v^Asn = n x 0.757 THz), 
implying that the SL is equivalent to a periodic SL with AlAs layers of thickness (Iaias and 
GaAs layers of thickness d^QaAs- Thus, the phonons can propagate through the random SL 
without seeing any of the structural disorder present in it. 

Narrow and wide band gaps coexist in the periodic double-layer SL made with the same 
unit blocks as the ones used for the random double-layer SL's. For these double-layer struc- 
tures, the location of the dips in the phonon transmission rate T (for the periodic SL) and 
those in (T) (for the random SL's) are somewhat different. This difference is much more 
pronounced when comparing T and < T > for the single-layer SL's. In the latter, the 
transmission dips in < T > do not appear at the frequencies within the narrow gaps present 
in the T of the periodic SL. We also note that the frequencies at which the perfect and 
almost perfect transmissions occur coincide with the v^ R > given above. 
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VI. LYAPUNOV EXPONENT 

Localized states can be characterized by the localization length £, or the Lyapunov exponent 
7, which is the reciprocal of the localization length, i.e., 7 = l/£. The Lyapunov exponent 
7 is defined by the average of the logarithmic decrement of the transmission rate as^l; 

7 = - lim -L<logT(L, W )>. (74) 

The Lyapunov exponent is subject to the following differential recursion relation^ for an 
arbitrary function f(z) of z(— l/T): 



!</> = <(2,-l)|> + W ,- l)g>. (75) 



The derivation of this relation is given in the Appendix. Putting f(z) = logz with the initial 
condition (f(z)) = 0ati = 0(L = 0or£^ 00), we find 

(-logT) = t. (76) 

From Eqs.(|74D and (76), the Lyapunov exponent 7 is given by 



7(w) =9T-TT- (77) 
2%l(^) 

Figure 3(a) shows the analytically obtained Lyapunov exponent 7 (oc I S (L — > 00)) versus 
frequency together with the numerical data for the single-layer random SL's. The numerical 
data are obtained for the random SL's with 1500 blocks for which the size effect is negligibly 
small. We see that 7 oscillates regularly and its peaks occur periodically at frequencies 
= n x 0.266 THz, which coincide with the Bragg frequencies of the periodic SL (the 
height of the peaks is the same, i.e., 1.158 x 10 _4 A _1 ). The Lyapunov exponent vanishes at 
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the resonance frequencies z/j n ( j = AlAs or GaAs) and varies parabolicly in the neighborhood 



of these frequencies 



7(f) 



8tc 2 R 2 d) 

d-AlAs + dcaAs Cj 




(78) 



In a finite system, phonons whose localization length is longer than the length of the 
system are regarded as extended. Thus, the critical localization length £ cr = 7" 1 = L (t cr = 
L/2) discriminates between localized and delocalized states. We can define the phonon pass 
band as a frequency region where the Lyapunov exponent 7 is smaller than the critical value 
7 cr . The width of the pass bands Au is estimated from Eq. (|78|) as 



Thus, the pass band width vanishes in an infinitely large SL as noted above. The average 
transmission at the band edges is readily estimated from Eq . (|67|) by putting t = t cr = 
L/£cr — 2. We find (T) = T cr = 0.26. Therefore, (T) is greater than 0.26 in the pass bands. 

Figure 3(b) shows the Lyapunov exponent 7 for the double- layer SL's, which exhibits 
quite a different frequency dependence from 7 in the single-layer SL's. In addition, the 
magnitude is significantly larger than 7 in the single-layer random SL's. These differences 
are attributed to the internal structure of the constituent blocks. In particular, the interfaces 
between any adjacent blocks are recognized by phonons, leading to a much shorter phonon 
mean-free-path as compared to single-layer SL's with the same number of unit blocks. This 
also leads to conspicuous dips in (T) in the double-layer SL's even when the number of blocks 
involved in a random SL is small (see Fig.2(b)). 

The parabolic behavior of the Lyapunov exponent 7 near the resonance frequencies and 
the size of a pass band width can be derived in a similar way as for the single-layer SL's but 
these calculations are not presented here. 




(79) 
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VII. UNIVERSALITY IN TRANSMISSION 



FLUCTUATIONS 



The phonons in a pass band are delocalized and their transmission rates are expected to fluc- 
tuate due to interference among multiply scattered phonons, just like the electronic conduc- 
tance fluctuations in mesoscopic transport^ 0. The standard deviation of the transmission 
AT, defined by 



can be readily calculated by using the formula derived from Eq . ([75]) , 

(T 2 ) = - i<T>. (81) 

Thus, the transmission fluctuation AT is also a function of t — L/£, and AT does not depend 
on the details of £, or the structure of the random SL's. This means that the relation between 
AT versus (T) should be universal, independent of the type of SL. 

Figure 4(a) shows AT versus (T) together with the numerical data for the single- and 
double-layer random SL's. Indeed, no difference can be seen in the numerical results of AT 
versus (T) and the analytical result coincides remarkably well with the numerical data. For 
localized phonons ((T) ~ 0), the fluctuations appear to be small and they increase with 
increasing (T). The maximum of AT is attained in the pass band at (T) ~ 0.4 where the 
mean-free-path becomes comparable to the system size, or L/t — 1.27 . 

Although the transmission fluctuations are small for localized phonons, the transmission 
rate cannot properly be averaged for these phonons. We plot the relative fluctuations of 
transmission in Fig. 4(b). The relative fluctuations diverge at (T) = and decrease with 
increasing (T). At the boundary between the localized and extended states, the relative 




(80) 
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fluctuations become almost unity and reduce to zero at (T) = 1. The relative fluctuations 
always exceed unity for the localized phonons, therefore, the transmission rate is not well- 
defined in the localized regime. 

It is known that, in a one-dimensional metal, the conductance or the resistance has the 
same properties as the phonon transmission rate discussed here. Anderson et al.0 applied a 
method of scaling transformation and found that the logarithm of the total resistance obeys 
a normal distribution for large sample lengths. 

In order to see the behavior of relative fluctuations of log T, we evaluate the ratio of the 
standard deviation of logT to (logT). With the use of Eq.([7|), the standard deviation of 
logT is given by 



Figure 5 shows the analytical and numerical results for the relative fluctuation of log T. Our 
analytical result indicates that the relative fluctuations become unity at (logT) = and 
decreases monotonically with increasing t. Thus, we conclude that logT, in contrast to T, 
is a well-defined quantity for localized phonons. 

The peculiar behavior of AT versus (T) can be explained in terms of the fluctuations 
of the Lyapunov exponent or the localization length plotted in Fig. 5. Here we note that 
(logT) is proportional to the Lyapunov exponent 7. The relative fluctuation of logT is, 
hence, equivalent to A7/7. For highly extended phonons with (T) ~ 1 where 7 is close to 
0, A7 has the same magnitude as 7 ( A7/ 7 ~ 1 from Fig. 5 ), i.e., A7 ~ 0. This is the 
reason why AT is small in the region (T) ~ 1. The relative fluctuation A7/7 decreases 
with increasing 7 (oc t). The relative fluctuation of the localization length, A£/£, is equal 
to A7/7 and hence, A£ is always smaller than £. In the extended region, the SL size L is 




(82) 
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shorter than £ but might become comparable to A£ because A£ ~ £. In this region, the 
transmission should be largely affected by fluctuations of the localization length or Lyapunov 
exponent. The fluctuation A£ is much shorter than the SL size L in the localized regime since 
A£ < ^ < L. Therefore, the fluctuation AT becomes quite small and vanishes at (T) ~ 0. 
However, we should keep in mind that the fluctuation A£ still affects the transmission in 
this regime because AT/ (T) diverges at (T) ~ as discussed above. 



In this section, we study the probability distribution of T in the random SL's. At the 
resonance frequencies, the transmission rate is unity and its distribution function should 
have a peak like 5(T — 1). From the analogy with the electronic localization problem, 
the probability distribution for localized phonons is expected to be Gaussian with respect 
to logT -1 . We are interested in the transition of the distribution with decreasing average 
transmission. We here introduce a probability density W(z,t) with which one can express 
the average of an arbitrary function f(z) of z = 1/T as 



Substituting Eq . (p3|) into ( ff5|) and integrating by parts, we obtain the Fokker-Planck equa- 
tion for W(z,t): 



VIII. 



LOG-NORMAL DISTRIBUTION OF 



TRANSMISSION 




(83) 



dW 



d_ 

dz 



[{2z - 1)W] + 



[z(z - 1)W] . 



(84) 



dt 



dz 2 



27 



Here the initial condition is W(z, t = 0) = S(z — 1). Equation ( p4|) was studied by Abrikosov 
and solved0 with the associated hypergeometric functions as 

W{z,t) = —/==] =exp [- [t/A + x 2 /t)\ dx, (85) 



where x = cosh -1 y/z. 

Figure 6(a) shows the statistical distribution of logT for the extended phonons (t = 0.47 
and 1.18, and the corresponding (T) are 0.67 and 0.42, respectively) and the critical states 
(t = 2, and correspondingly (T) = T cr = 0.26) together with the numerical data. The 
similar plot for the localized phonons (t = 8.22, and correspondingly (T) = 0.02) is given in 
Fig. 6(b). The ordinate has been taken as W(z,t) z because 

dz 

W(z,t)dz = W(z,t) — — - dlogT = -W(z,t) z d(logT). (86) 
ctlogT 

The distribution for the extended phonons has a peak at the origin logT = 0, suggesting 
a Poisson distribution. The peak height is lowered with decreasing (T). Simultaneously, 
the FWHM increases and the shape of distribution curves becomes convex upward, but the 
position of the peak stays unchanged (logT = 0). For the strongly localized states, i.e., 
(T) T cr , a Gaussian distribution is attained. The distribution W(z,t) z for the localized 
states is analytically derived from Eq.(RH) as0 



W(z,t)z ~ — ^=exp 
2v nt 



(logT + tf 



, t»l, (87) 



4f 

and the peak is located at a finite value of logT, i.e., logT = —t. The solid curve of Fig.6(b) 
for t = 8.22 agrees well with Eq.(0)- The FWHM obtained from Eq.(g7D is \/2t. From 
these quantities, the relative fluctuations of logT for the strongly localized states can be 
calculated again as \/2t/t = \j2/t, which coincides with the asymptotic formula (valid at 
t ^> 1) derived analytically from Eq.(^) (see Fig.5). 
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To conclude, the transmission rate in the localized regime obeys the log-normal distri- 
bution and a Poisson distribution applies to highly extended phonon states. 

IX. SUMMARY AND CONCLUSION 

Based on the Green's function method, we have studied the transmission of phonons prop- 
agating along the growth direction of random SL's. 

Taking account of the phase shift of phonons associated with the forward scattering due 
to the mass density fluctuations, we have derived an explicit expression for the average 
transmission (T) as a function of the scaling parameter t = L/£, i.e., the size of the SL 
divided by the mean-free-path of phonons limited by backscattering. The mean-free-path £ 
reflects the structural properties of the SL and is expressed by the average of the squared 
SL structure factor, (\S S l(L, uj)\ 2 ). We can readily evaluate the average transmission (T) for 
a random SL, by calculating analytically the mean-free-path £. We applied the analytical 
expression for the average transmission (T) to two types of SL's with different unit-block 
structures, i.e., single- and double-layer random SL's. The average transmission rates show 
dips and peaks due to coherent backscattering and resonances of phonons depending on the 
structure of the SL. The coincidence of the analytical results with the numerical data based 
on the transfer matrix method verifies the accuracy of our theory. The universality relation 
of AT versus (T) has also been established explicitly for both types of SL's. Finally, we 
have examined the statistical distribution of the phonon transmission and have confirmed a 
log-normal distribution for the localized phonons. 

The most important result we have found in the present work is that the phonon-mean- 
free-path (and eventually (T)) can be calculated from (\Ssl(L,u)\ 2 ). This enables us to 
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study quantitatively the average transmission for a given set of SL's with simple calculations 
of the mean-free-path based on the prescribed statistical rules for the building blocks. 

The existence of resonances, in the SL system we have considered, yields an interest- 
ing structure in the phonon transmission spectrum, such as the alternating regions, in the 
frequency domain, of localized and extended states. By modulating the thickness of the 
constituent layers, we can readily produce resonances in the sub-THz frequency region, which 
are accessible experimentally. In pure semiconducting materials, the phonon mean free path 
limited by the bulk elastic scattering caused by mass defects (foreign and isotopic impurities 
etc.) is 0.1 to 1 cm, or much longer at frequencies below ITHz. Also, inelastic scattering 
of phonons is much weaker at low temperatures. Hence, semiconductor hetero-structure 
random SL's would permit the observation of the effects of high-frequency phonon multiple- 
scattering induced by the mass-density fluctuations in the layered structures. 

Through this work we have assumed the same stiffness constants for the SL constituent 
materials. However, this does not limit the applicability of our results because the fluctua- 
tions in the stiffness constants are effectively incorporated into the mass-density fluctuations 
in the present formulation. Thus, our results are readily applicable to SL's with various 
combinations of materials, and also to the study of propagation of third-sound waves in Hell 
films adsorbed on disordered substrates^. 
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APPENDIX: 
Derivation of Differential Recursion Relation 

Consider an arbitrary function f(z) expanded as 

f(z) = J2a n z n , (Al) 

n 

where z = l/T(L,u) = [S(L, 0)] u [S(L, 0)] 22 . Substituting Eq. (pl|) into (A.l) and expanding 
the function with respect to n and n*, an ensemble average of a term z n is given up to 0(|k| 2 ) 
by 



(*T> = (4) + (\^}{n(z^(2z - 1)) + n(n - l)(z^\z - 1))}, (A2 



where z\ = z(L) and Zq = z(L — A). Taking the limit of A — > and dividing (A. 2) by (|k| 2 ), 
we obtain the following differential equation for the (z n ) term 

_(,..) = ( _ (2z _ !), + <_,(,_!». (A3 ) 
Thus, the differential recursion relation for an arbitrary function f(z) is derived as 

|c/W)-<^P.-D> + (^(.-i)). <M) 

This equation was originally derived by AbrikosovEi for the study of the electrical conduc- 
tivity in disordered systems. 
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Figure captions 



Fig. 1 Schematic representations of the building blocks for (a) single-layer SL's and (b) 
double-layer SL's made of two kinds of materials A and B. In the numerical calculations 
AlAs and GaAs are assumed for the materials A and B, respectively. 

Fig. 2 Average phonon transmission rate, (T), versus frequency, v, for (a) single-layer 
random SL's with 1000 building blocks (N = 1000) and (b) double-layer random SL's with 
100 building blocks (N = 199). The solid lines show the analytical results and the open 
circles are the numerical data calculated from the transfer matrix method. Each data point 
is obtained by averaging over 100 realizations of SL randomness. The parameters used are: 
d/aAs — dcaAs = 34A for the single-layer random SL's used in (a); and d^AiAs — d 2 ,AiAs = 
dAiAs = 17 A, d\ : G a As = 20A and d2,GaAs = 42A for the double-layer random SL's used 
in (b). The mass densities are Paias — 3.76 g cm~ 3 and pcaAs = 5.36 g cm~ 3 ; and p = 
5.915 x 10 11 dyn cm~~ 2 . The slight shifts of resonance frequencies {vg}as,u an< ^ u MAs,n) f rom 
those in Ref.3 are attributed to the average stiffness constant p used in the present model. 
The dotted line shows the transmission rate of a periodic SL with the same constituent 
blocks. 

Fig. 3 Lyapunov exponent 7 (the reciprocal of the localization length £) versus fre- 
quency, v, for (a) the single-layer random SL's and (b) the double-layer random SL's. The 
solid lines show the analytical results and the open circles are the numerical data. The 
numerical results are obtained for (a) the single-layer random SL's with iV = 1500 and (b) 
the double-layer random SL's with iV = 199, respectively. The data are averaged over 100 
possible realizations of SL randomness. 
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Fig. 4 (a) Standard deviation of the phonon transmission rate, AT, versus (T). The 
continuous line denotes the analytical result and the open squares and circles are the numer- 
ical data obtained for (a) single- (open squares) and (b) double- (open circles) layer random 
SL's, respectively. Each data point is obtained by averaging over 100 different SL's. (b) 
The relative transmission fluctuation AT/(T) versus the average, (T). The solid line is the 
analytical result and open circles are numerical data for the two types of random SL's con- 
sidered. T cr = 0.26 (L = 21) indicates the boundary between the extended and localized 
regimes for phonons. 

Fig. 5 Relative fluctuation of logT versus — (logT) = t. The solid lines show the 
analytical results and the open circles are the numerical data for the two types of random 
SL considered. The dashed line is the asymptotic formula (t/2)~ 1 / 2 derived analytically for 
t > 1. 

Fig. 6 Probability distribution of W(z,t) z versus —logT togetherwith the numerical 
data for (a) extended phonons ((T) > T cr = 0.26) with (T) = 0.67 (solid line and circles), 
0.42 (dashed line and squares) and critical states with (T) = T cr , i.e., L = £ cr (dot-dashed 
line and triangles), and for (b) localized phonons with (T) = 0.02 (solid line and squares). 
The arrows indicate the positions of t = (—logT). Note that z = 1/T. The numerical data 
are obtained for an ensemble of 1000 single-layer SL's by calculating T at a single frequency 
giving the chosen value of (T). 
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